Determination of transverse momentum depen- 
dent gluon density from HERA structure func- 
tion measurements 



Hannes Jung 1 ' 2 , Francesco Hautmann 3 

1 DESY, Notkestrafie 85, 22607 Hamburg, Germany 
2 CERN, 1211 Geneve 23, Switzerland 

3 Theoretical Physics Department, University of Oxford, Oxford 0X1 3NP, GB 
DOI: will be assigned 

The transverse momentum dependent gluon density obtained with CCFM evolution is 
determined from a fit to the latest combined HERA structure function measurements. 



1 Introduction 



The combined measurements of the structure function at HERA 1 allow the determination 
of parton distribution functions to be carried out to high precision. While these data have 
been used to determine the collinear parton densities, the transverse momentum distributions 
(TMD) or unintegrated gluon distributions were only based on older and much less precise 
measurements [2J|3]- 

In high energy factorization 0] the cross section is written as a convolution of the partonic 
cross section a(k t ) which depends on the transverse momentum k t of the incoming parton with 
the /c t -dependent parton density function A (x, kt,p): 

—d 2 k t a(-,k t )A(x,k t ,p) (1) 

z z 

where p is the factorization scale. The evolution of A{x 1 k t ,p) 
can proceed via the BFKL, DGLAP or via the CCFM evolution 
equations. Here, an extension of the CCFM evolution is applied 
(to be also used in the parton shower Monte Carlo event generator 
CASCADE [5]) which includes the use of two loop a s as well as 
applying a consistency constraint [SJ H] in the g — > gg splitting 
function [9]: 

'(1-z) , z(l-z) 




P gg (z,p,k t ) = a s (k 2 ) 



A ns + a s (p ) 
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Figure 1: Gluon branching 
z z(l — z) " 
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with A„ s being the non-Sudakov form factor. The consistency constraint is given by [5] (see 
Fig.0: 

q 2 t < V >-L ( 3 ) 
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2 Evolution 



Since the CCFM evolution cannot be easily written in an analytic closed form, a Monte Carlo 
method, based on [TOMTT], is used. However, the Monte Carlo solution is time consuming, and 
cannot be used in a straightforward way in a fit program. For a realistic solution, first a kernel 
A (x" , k t ,p) is determined from the MC solution of the CCFM evolution equation, and then is 
folded with the non-perturbative starting distribution Aq(x): 



xA(x,k t ,p) = x j dx' j dx"A (x)A(x",k u p)6(x' -x" - x) (4) 

dx' [ dx"A (x)A (x", k uP ) -8{x" - ~) (5) 

J X X 

o'A (x')-^A^,k u p) (6) 



dx'. 

The kernel A includes all the dynamics of the evolution, Sudakov form factors and splitting 
functions and is determined in a grid of 50 ® 50 ® 50 bins in x,kt,p. 

The calculation of the cross section according to eq. ([I]) involves a multidimensional Monte 
Carlo integration which is time consuming and suffers from numerical fluctuations, and cannot 
be used directly in a fit procedure involving the calculation of numerical derivates in the search 
for the minimum. Instead the following procedure is applied: 



a r (x,Q ) = / dx g A(x g ,k t ,p)a(x,x g ,Q ) (7) 

dx g dx' dx"A (x')A(x", k u p)- &(x, x g , Q 2 ) ■ S(x' x" - x g ) (8) 

dx' dx"A {x') ■ A(x", k t ,p)- a(x, x' x" , Q 2 ) (9) 

r 1 

dx'Ao(x')- / dx"A{x'\k tl p)-a{x,x'x",Q 2 ) (10) 

J x I x' 

dx'A (x')-a(x/x',Q 2 ) (11) 

J X 

Here, first a(x',Q 2 ) is calculated numerically with a Monte Carlo integration on a grid in x 



for the values of Q 2 used in the fit. Then the last step (i.e. eq.(lll) is performed with a fast 
numerical gauss integration, which can be used in standard fit procedures. 

The fit to the HERA structure function measurements is performed applying the heraf itter 
program [I] [T2 H2] to determine the parameters of the starting distribution _4o at the starting 
scale Qq: 

xA (x, k t ) = Nx- B * ■ (1 - x) c ° (1 - D g x) (12) 

3 Fit to HERA structure function 



The parameters N,B g ,C g ,D g in eq.(12) are determined from a fit to the combined structure 



function measurement pQ in the range x < 0.01 and Q 2 > 5 GeV. In addition to the gluon 



2 



DIS 2012 



induced process 7*5* — > the contribution from valence quarks is included via 7*17 — > g 
using a CCFM evolution of valence quarks as described in [2]. The results presented here 
are obtained with the heraf itter package, treating the correlated systematic uncertainties 
separately from the uncorrelated statistical and systematic uncertainties. To obtain a reasonable 
fit to the structure function data, the starting scale Qq as well as A qcc i has been varied. An 
acceptable x 2 / n df could only be achieved when applying the consistency constraint cq.([3|: 
without consistency constraint the best x 2 jndf ~ 14 — 28, depending on which form of the 
splitting function is used. With consistency constraint and the splitting function eq.([2]) the 
best fit gives x 2 / n df ~ 1.5 for Q — 1.8 GeV and A qc d — 0.17 GeV at 71/ = 4 flavours. It 
has been checked, that the \ 2 l n df does not change significantly when using 3 instead of 4 
parameters for the initial starting distribution Aq. 

In fig{2] the resulting unintegrated gluon density JH-setO is shown for 2 values of p 2 com- 
pared to set AO [T5] , 



p = 0.25E+O2 GeV" 
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Figure 2: Unintegrated gluon density JH-setO for p 2 = 25 GeV 2 (left) and p 2 = 10 5 GeV 2 
(right) as a function of x for different values of k 2 and as a function of k 2 for different values of 
x compared to set AO [T5] 

The uncertainties of the pdf are obtained within the heraf itter package from a variation 
of the individual parameter uncertainties following the procedure described in |16| applying 
A\ 2 = 1. The uncertainties on the gluon are small (much smaller than obtained in standard 
fits), since only the gluon density is fitted. The uncertainty bands for the gluon density are 
shown in fig. |3jleft). In fig. fright) the prediction for 6-jet cross section as calculated from 
Cascade [5] using the gluon density described here (labeled as setO) is shown together with a 
prediction using an older set (labeled as set A |15j ) in comparison with a measurement from 

cms [13. 
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Figure 3: (left): Uncertainties of the uPDF at p 2 = 25 GeV 2 . (right): Cross section of 6-jet 
production as a function of pt for different bins in y as measured by CMS |17j compared to 
predictions from Cascade [5] using the unintegrated gluon density described here 



the meeting. 
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